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ABSTRACT. An ice-surface temperature retrieval algorithm for the Greenland 
ice sheet was developed using NOAA 1 1 thermal radiances from channels 4 and 5. 
Temperature, pressure and humidity profiles, cloud observations and skin 
temperatures from the Swiss Federal Institute of Technology (ETH) camp, located 
at the equilibrium line altitude at 49° 17' W, 69°34'N, were used in the LOWTRAN 
7 model. Through a statistical analysis of daily clear sky profiles, the coefficients that 
correct for the atmospheric effects were determined for the ETH-Camp field season 
(May to August). Surface temperatures retrieved by this method were then 
compared against the in situ observations with a maximum difference of 0.6 K. 
The NOAA 1 1 narrow-band planetary albedo values for channels 1 and 2 were 
calculated using pre-launch calibration coefficients. Scattering and absorption by the 
atmosphere were modelled with LOWTRAN 7. Then, narrow-band albedo values 
for the AVHRR visible and near infrared channels were compared with in situ high 
resolution spectral reflectance measurements. In the visible band (580-680 nm), 
AVHRR-derived narrow-band albedo and the in situ measurements corrected with 
radiative transfer model LOWTRAN 7 showed a difference of less than 2%. For the 
near infrared channel (725-1 100 nm) the difference between the measured and 
modelled narrow-band albedo was 14%. These discrepancies could be either the 
result of inaccurate aerosol scattering modelling (lack of the in situ observation), or 
the result of sensor drift due to degradation. 


INTRODUCTION 

Past and ongoing research indicated that several polar 
surface properties can be mapped through a combination 
of multispectral satellite data. These properties include 
spectrally integrated surface albedo, surface temperature 
and short- and long-wave radiation balance, which in 
turn is an important measure for the discussion of the 
long-term climate stability. The Greenland ice sheet 
provides an ideal case for the calibration of satellite- 
derived parameters due to its large size and different 
homogenous snow and ice surfaces over large areas. 
Further, the study of seasonal and interannual variations 
of the above mentioned parameters provide an ideal data 
set to determine the “noise level” for future climate trend 
studies. Here the technique and some preliminary results 
are presented to estimate ice surface temperatures and 
surface albedo values from different spectral channels of 
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DATA AND METHOD 

AVHRR calibration and navigation 

The advanced very high resolution radiometer 
(AVHRR) data used in this study were obtained from 
National Oceanographic and Atmospheric Research 
Laboratories (NOARL). This is the local area coverage 
(LAC) data in level lb format from the NOAA-11 
satellite, and includes 5 channels: 2 visible and 3 infrared. 
The level lb format includes information for the 
navigation and calibration of the images (Lauritson and 
others, 1979), but this information was re-done on site to 
gain more control over the processing of the data. The 
navigation of the AVHRR data can be used to better 
depict the actual motion of the satellite (Baldwin and 
Emery, 1993). A small registration error is present in the 
resulting images that are removed by linearly nudging the 
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are then compared with each other. Registration errors 
can then be reduced to pixel. The navigation of the 
data was mapped to a polar stereographic projection with 
a 1 km per pixel resolution. The resolution of the scan spot 
from the AVHRR instrument varies from 1.1 km at nadir 
to about 6.8 km at a scan angle of 55°. This is not a linear 
function; a scan angle of 30° still results in a pixel 
resolution of about 1.5 km. 

Calibration of the AVHRR imagery is channel 
dependent. Calibration of the reflected solar energy 
“visible” channels, 1 and 2, is done by using pre-launch 
values particular to each satellite. These values (Kidwell, 
1991) convert the instrument counts into per cent albedo. 
Solar zenith angles are computed during the navigation 
procedure for each pixel in the image and corrections 
made to the per cent albedo. Channel 3 is a combination 
of reflected and emitted energy during daylight condit- 
ions and was not used in this paper. The emitted infrared 
channels, 4 and 5, require several steps for calibration. 
First, continuous on-board sensing of space and black- 
body targets every few scan lines provides information for 
computing slope and intercept values for the infrared 
channels (Lauritson and others, 1979). Energy values are 
calculated, and then using Planck’s equation with pre- 
launch wavenumbers, brightness temperatures can be 
found. Non-linear corrections that are dependent upon 
the scene temperature are then applied. These non-linear 
corrections were on the order of -1.0° for the observed 
area. 

Ice-surface temperature 

A common approach for estimating surface temperatures 
is to relate satellite data to surface temperature observa- 
tions with a regression model. The radiative transfer 
equation can be applied to model satellite radiances and 
brightness temperatures. A common approach to retrieve 
surface temperatures is to use the split-window infrared 
AVHRR channels 4 and 5 (10.5-1 1.5 /im, 1 1 .5-12.5 /im) 
as was applied for open-ocean areas (McClain and others, 
1985; Minnett, 1990) as well as sea-ice areas (Key and 
Haefliger, 1992). To correct for atmospheric attenuation 
of the satellite measured brightness temperatures over the 
Greenland ice sheet the surface temperatures are 
regressed against modelled AVHRR brightness temper- 
atures as described by Key and Haefliger (1992). Atmos- 
pheric attenuation was derived with the radiative transfer 
model LOWTRAN 7 (Kneizys and others, 1988) using 
radiosonde profile measurements from the ETH-Camp 
(Federal Institute of Technology (ETH), Department of 
Geography (Fig. 1)) as input data. A total of 78 
radiosondes were launched from 8 June through 29 
August 1990, and 194 from 13 May through 21 August 
1991. The profiles include pressure, temperature and 
relative humidity for every second up to a height of 20- 
25 km (Fig. 2). Continous incoming and outgoing 
longwave radiation measurements and daily cloud 
observations (7 times each day) were also available from 
the ETH Camp. The longwave radiation at the ETH 
Camp was calculated as the difference between the all- 
wave radiation, measured by pyrradiometers, and the 
short-wave radiation, measured by pyranometers, and 



Fig. 1. Overview of Disco Bay in West Greenland 
(560 x 560km). The rectangular area around the ETH 
Camp ( 152 x 227 km ) is the test area to which the 
coefficients are applied ( dotted lines — coast , solid line = 
ice margin). 


ETH-Camp 90-6-27 23:52 



Temperature (C) / Windspeed (m/s) / Rel. HumWity {%) 


Fig. 2. Radiosonde profile measurements of temperature , 
relative humidity , and wind speed from the ETH Camp , 
located at the equilibrium line altitude on the Greenland ice 
sheet (49° 17 69°34'N). 


The pyrradiometers were calibrated over a melting ice 
surface. Only situations with less than | cloud cover were 
used in the analysis, which are referred to as “clear sky” 
conditions. A total of 90 profiles remained with clear sky 
conditions. 

The surface temperature at the camp is calculated 
from the outgoing longwave irradiance with the Stefan- 
Boltzmann law: R = <tT 4 , assuming the snow surface is a 
black-body. In general the longwave irradiance measur- 
ements have an error of ±10Wm~ 2 (±3%) or even 
better in case of the Greenland measurements (personal 
communication from T. Konzelmann). To every meas- 
ured outgoing longwave irradiance 3% are added and 
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created out of each of the 90 original clear sky ones, where 
the surface temperatures are the minimum 



and the maximum 



calculated temperatures, respectively. This results in a 
sample size of 270 clear sky profiles. 

To eliminate the influence of the atmosphere on the 
signal reaching the satellite, the measurements of at least 
two spectral bands are required. A good overview of the 
theoretical algorithms for satellite-derived SSTs is 
presented in Barton and others (1989). The method 
used in this study is the same as described in Key and 
Haefliger (1992) and reads as following: 

T ice - a + 6 T 4 + cT 5 + d[(T A - T 5 )s(0)]. (1) 

T 4 and T 5 are the satellite-measured brightness temper- 
atures in Kelvin and 0 is the sensor scan angle (0-55°). 
The coefficients are determined through a least-squares 
regression procedure. Surface temperatures are regressed 
against modelled brightness temperatures. LOWTRAN 7 
is used to simulate radiances in the AVHRR thermal 
channels 4 and 5, which are then converted to brightness 
temperature. 

The coefficients for the estimation of ice surface 
temperature (1ST) according to Equation (1) are shown 
in Table 1. The root mean square (RMS) error of 0.3 K is 
well within the error of the outgoing longwave measure- 
ments which is + 1 K for the Greenland measurements. 
Table 2 shows the calculated 1ST for the AVHRR pixel 
representing the location of the ETH Camp. Due to the 
uncertainty of the exact location (AVHRR geolocation 
error), the temperature of 3 x 3 pixels are calculated 
and compared to the in situ measurement of the camp. 
The comparison of AVHRR-derived ice surface tempera- 
tures with the in situ measurements show differences of up 
to 0.6 K. The regional variation of ice surface tempera- 
tures around the ETH Camp are shown in Figure 3 (see 
also Figure 1 for overview of area). 

Planetary and surface narrow-band albedo 

In the presented case study, the planetary narrow-band 
albedo for AVHRR channels 1 and 2, measured on 23 
May 1991, will be compared against ground-based values 
measured with a portable spectrometer. The Geophysical 
Research Corporation (GER) portable field spectrometer 
measures reflectance values for 875 different spectral 
bands between 300 and 2500 nm (Table 3). The spectro- 
meter was regulated with a laptop computer (SNAP 
1 + 1) which was also used for data acquisition. The field- 
of-view of the optical head is 7° in horizontal and 1.6° in 
vertical direction (rectangular area of sensitivity given by 
the photodetector array). An integration sphere was 
mounted in front of the optical head for hemispheric 
measurements which could be directed towards the zenith 


Table L Coefficients and RMS error based on NO A A 11 
A VHRR channels 4 and 5 


a b c d RMS 

-4.257151 3.473293 -2.470502 —0.141503 0.320 


Table 2. Measured and calculated AVHRR ice-surface 
temperature (1ST) at the ETH Camp 


Date 

1990 

Measured 1ST 
K 

Calculated 1ST 
K 

Delta 1ST 
K 

21 June 

271.1 

271.6 

+ 0.5 

23 June 

271.3 

271.0 

-0.3 

24 June 

271.5 

271.1 

-0.4 

27 June 

271.3 

270.8 

-0.5 

28 June 

271.3 

271.1 

- 0.2 

29 June 

270.2 

270.7 

+ 0.5 

1 July 

271.9 

271.3 

- 0.6 

2 July 

271.1 

271.5 

+ 0.4 


Table 3. Specification for the GER- spectrometer 


Range 

nm 

Resolution 

nm 

Number of Channels 

300- 600 

1.5 

199 

600-1050 

1.7 

266 

1050-1800 

3.0 

245 

1800-2500 

4.2 

165 


measure the spectral reflectance of the snow surface. 
The spectral reflectance of the snow surface is shown in 
Figure 4. The snow was fresh, less than one day old, with 
grain sizes between 0. 1-0.2 mm in diameter. The 
prevailing air temperature was well below freezing. In 
addition to the spectrometer measurements broadband 
global radiation, and diffuse and reflected shortwave 
radiation were measured with a set of pyranometers as 
part of the ongoing radiation balance experiment at the 
ETH Camp. The reflected shortwave radiation was 
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< 280.0 

277.5 

276.5 

275.5 

274.5 

273.5 

272.5 

271.5 

270.5 

269.5 

268.5 

267.5 

266.5 

265.5 

264.5 

263.5 

> 260.0 


p ig. 3. Ice-surface temperatures around the ETH camp (see Fig. 2). The variation from day to day is shown. The image 
of 26 June contains an error in a few scan lines } visible in the lower right comer. 


broadband albedo values derived from the pyranometers 
were 84.3% for the lower instrument, and 82.2% for the 
instrument at 27 m. The slightly lower albedo value for 
the later can be explained by the larger field-of-view 
which included part of the field camp as well as surface 
roughness in the scale order of meters and larger. The 
spectrally-integrated reflectance derived from the spectro- 
meter is 82.9% which shows the level of accuracy 
(6 — 0.7%) which can be expected using two different 
kind of instruments. 

The narrow-band albedo measured by the AVHRR 

( • 1 ■ .i r. . « . 


input value C and the pre-launch calibration coefficients 
as follows (Kidwell, 1991): 

A = S(i)C + I(i ) (2) 

with 5(1) = 0.095, 5(2) = 0.1061, 7(1) = -3.8, 1(2) = 
-3.6 for NOAA 11. 

This calibration takes the spectral response of the 
channels (t = 1,2) into account by multiplying the 
spectral irradiance viewed at the top of the atmosphere 
with the spectral response function. NOAA 1 1 satellite 
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Greenlond Station: New Snow, Groin Size - 0.1 -0.2 mm (5/23/91) 



Wavelengths (nm) 

Fig . 4. Spectral reflectance of new snow (grain size 0.1- 
0.2 mm ) , pre-melt period. 


1991, 15.11 GMT with a satellite zenith angle relative to 
the camp of 27. 37° and a solar zenith angle of49.16°. The 
narrow-band planetary albedo values for AVHRR 
channels 1 and 2 are listed in Table 4. The accuracy of 
the AVHRR-derived narrow-band albedo values are 
discussed in the following. 

The atmospheric absorption and scattering was 
modelled using the radiative transfer code LOWTRAN 
7. Two model runs were performed: (a) using the in situ 
radiosonde profile measurements of temperature and 
humidity as input data, and (b) using the standard 
subarctic atmosphere. Figure 5 depicts the solar irrad- 
iance at the top of the atmosphere and the irradiance as 
modelled with LOWTRAN 7 at the surface using 
radiosonde profile measurements. The narrow-band 
albedo at the top of the atmosphere a p was then 
calculated by integrating over the AVHRR spectral 
band-width (Ai,A 2 ), as following: 

/[5(A) 0j (A)r(A)a g (A)]dA , . 

“ P /[5(A) 0, (A)]dA 1 ' 

where S(A)is the solar irradiance at the top of the 
atmosphere, 0i(A)is the spectral response function of the 
AVHRR channel (i = 1,2), a g ( A) is the narrow-band 
reflectance at the ground, and r(A) is the calculated 
spectral transmission for two path lengths: (a) between 
the top of the atmosphere and the surface, and (b) 
between the surface and the satellite. The integration was 
done for both AVHRR channels with (a) radiosonde 
profile data and (b) standard sub-Arctic atmosphere 
(Table 4). The difference between measured and 
calculated narrow-band planetary albedo for channel 1 
is 0.5% using radiosonde profile data in the LOWTRAN 
7 model. For the standard sub-Arctic atmosphere, the 
difference increased to 2.9%. However, for AVHRR 
channel 2 the calculated and measured narrow-band 
albedo values differ by approximately 9%. Surface 
narrow-band albedo a g can also be derived from 


Table 4. Narrow-band albedo values for AVHRR 
channels 1 and 2 as measured by A VHRR at the top of 
the atmosphere (dp = planetary narrow-band albedo ) } 
and calculated using ground-based spectrometer values a g 
and radiative transfer modelling 


AVHRR a p Calculated a p Calculated a p a g 

Channel with radiosonde with standard 

atmosphere 

% % % % 


1 (580-680 nm) 68.8 67.5 65.5 91.6 

2 (725—1 100 nm) 56.1 65.2 64.5 81.0 


AVHRR 1 AVHRR2 



Fig. 5. Solar irradiance outside the atmosphere , and at the 
surface as modelled by LOWTRAN 7 radiative transfer 
model with the in situ radiosonde data for 23 May , 15.11 
GMT. 


Table 5. Narrow-band surface albedo values a g compared 
against calculated surface narrow-band albedo from 
A VHRR channels 1 and 2 using transmission coefficients 
n ( sun-surface ) y and 7*2 (surface-satellite) 


AVHRR 

a g 

Calculated 

Calculated 

n 

r 2 

Channel 

a g with 

a g with 





radiosonde 

standard 

atmosphere 




% 

% 

% 



1 (580-680 nm) 

91.6 

93.3 

95.1 

0.825 

0.895 

2 (725-1 100 nm) 

81.0 

69.3 

70.3 

0.878 

0.922 
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[^(suii -surface) (*^) ^(surface -satellite) (^)J 


The transmission values were calculated with the 
radiative transfer model which takes the intrinsic atmos- 
pheric reflectance due to the Rayleigh and aerosol 
scattering of the standard sub-Arctic atmosphere into 
account. The surface narrow-band albedo values derived 
from the AVHRR planetary albedo are given in Table 5. 
The differences have the same order of magnitude for the 
calculated planetary narrow-band albedo which is not 
surprising as both computations use the LOWTRAN 7 
model with identical input parameters. 

So far, we have neglected the bidirectional reflectance 
(BR) of a snow surface as a possible error due to large 
scanner angle and large solar zenith angle (Steffen, 1987). 
It can be assumed that the snow has an isotropic 
reflectance for the small satellite zenith angle encoun- 
tered during this study. For satellite zenith angles of 50° 
or larger, the BR of snow has to be taken into account, 
especially for wet snow surface. The in situ measurements 
of aerosol distribution with height were not available. 
This might be an additional error source. Aerosol content 
can vary with time and space significantly which can 
reduce the reflectance by several per cents. If we assume 
that the radiative transfer model is correct, we can 
speculate that a major sensor drift of AVHRR channel 2 
has occurred since launch of NOAA 1 1 on 8 November 
1988. For the visible channel of NOAA 9 degrading rates 
of up to 6% per year were found (Staylor, 1990). This 
case study should be judged as a preliminary analysis, as 
more cases have to be compared to get a statistically 
significant error analysis. Further, we have to find a 
means to validate the LOWTRAN 7 derived spectral 
transmission values, in order to separate the modelling 
error from the error caused by the sensor drift. 


SUMMARY AND CONCLUSION 

The relationship between clear sky NOAA 1 1 thermal 
radiances and the measured surface temperature of the 
Greenland ice sheet is examined through forward 
calculations of radiative transfer equation. Pressure, 
temperature and relative humidity profiles and cloud 
observation from the ETH Camp are used. Using the split 
window channels 4 and 5 and scan angle, the RMS error 
in the estimated 1ST is 0.3 K. 

Narrow-band surface albedo values can be estimated 
from AVHRR channel 1 and 2 measurements with an 
accuracy of 2% and 14% respectively, using LOWTRAN 
7 radiative transfer modelling and radiosonde profile 
information. More case studies are needed to discuss this 
technique in more depth and to derive statistically 
significant values. 
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